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Abstract 


Divergence  and  vorticity  are  computed  for  a  small  area  of  about  l/2 
square  mile  from  the  changes  in  area  and  orientation  of  a  triangle  formed 
by  three  pilot  balloons  as  they  ascend.  Mean  local,  seasonal  ascention 
rates  were  previously  determined  at  0300  Z  and  1^00  Z  for  the  purpose  of 
making  these  divergence  and  vorticity  computations.  The  magnitudes  and 
vertical  distributions  of  the  divergence  and  vorticity  are  discussed. 
Levels  of  least  divergence  are  located  at  about  6000  ft  and  13,000  ft, 
showing  close  agreement  with  an  earlier  study  which  dealt  with  a  scale 
of  motion  compatible  with  the  normal  upper  wind  network.  In  one  triple¬ 
balloon  run  there  is  strong  evidence  that  a  narrow  zone  (no  more  than 
3000  ft  wide)  may  exist  in  a  horizontal  plane  across  which  the  wind  direc¬ 
tion  changes  about  130*. 


1.  Introduction 


In  Part  1  of  this  study  [1]  mean  local  ascentional  rates  were  deter¬ 
mined  for  the  thirty-gram  balloon  at  the  times  0300  Z  (night)  and  1500  Z 
(day)  at  Tallahassee,  Florida  in  the  summer  of  195U.  The  details  are  de¬ 
scribed  in  [1]*  suffice  it  to  say  here  that  both  the  day  and  night  mean 
rates  were  found  to  be  significantly  greater  than  the  normally  assumed 
rate  after  the  first  5  or  i  minutes  of  ascent  (see  fig.  2  in  [1]).  The 
mean  local  rates  were  greater  than  the  assumed  rate  by  as  much  as  9  per 
cent  and  the  systematic  nature  of  the  difference  caused  the  height  of  the 
balloon  at  higher  levels  to  be  considerably  different  from  that  normally 
assumed.  For  example,  at  the  end  of  thirty  minutes  the  balloon  would  be 
found  at  18,012  ft,  18,602  ft  and  19,061;  ft  using  the  assuned  rate,  the 
mean  local  day  rate  and  the  mean  local  night  rate,  respectively. 

For  the  second  phase  of  the  study  it  is  important  to  know  as  accu¬ 
rately  as  possible  the  height  of  the  balloons  at  any  given  time  and  the 
rate  of  ascent  and  changes  in  this  rate.  This  was  the  reason  for  deter¬ 
mining  mean  local  day  and  night  rates  and  these  mean  rates  are  used  in 
determining  the  divergence  and  vorticity. 

2,  The  method 

A  triangle  of  theodolite  stations  was  set  up  (see  [1])  with  the  legs 
of  the  triangle  being  roughly  one  mile  in  length.  Three  balloons  released 
simultaneously  (one  at  each  site)  describe  a  triangular  area  at  each  min¬ 
ute  reading.  The  changes  in  size  and  orientation  of  this  area  can  then 
be  expressed  as  the  divergence  and  vorticity  of  the  flow  in  this  small 
region. 

For  the  horizontal  velocity  divergence  we  can  write. 


(1) 
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where  x  and  y  are  the  sides  of  area  A,  t  represents  time  and  u  and  v  are 
the  wind  speed  components  in  the  east  (x)  and  north  (y)  directions.  At 
time  zero,  A  is  the  area  of  the  fixed  triangle  of  theodolite  stations. 

The  area  at  the  end  of  each  minute  interval  can  be  obtained  from 
the  determinant 

xi  yi 1 

At  -  1/2  x2  y2  1  (2) 

x3  y3  1 

Fig.  1  shows  the  release  sites 

/ 

1,  2,  and  3  and  the  orientation 
of  the  coordinate  system;  also 
the  initial  area  in  ft^  and  the 
coordinates  of  the  sites  in  ft. 

The  slight  difference  in  height 
of  the  three  sites  is  neglected 
and  the  initial  triangle  is 
assumed  flat  for  the  divergence 
and  vorticity  computations.  The 
error  caused  by  this  assumption 

Fig.  1.  Coordinates  of  the  three 

is  of  the  order  of  a  fraction  of  balloons  at  release  time. 


the  horizontal  range  to  a  point  on  the  surface  directly  below  the  balloon. 
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is  given  by  *  zn  cot  en  in  which  z  is  the  height  of  the  balloon  above 
the  surface  and  e  is  the  elevation  an^lej  and  where  the  subscript  t  is 
the  time  of  reading  after  release,  in  minutes,  the  subscript  zero  refers 
to  the  release  tine,  and  the  subscrijh.  n  is  equal  to  1,  2,  or  3,  depending 
upon  which  release  site  is  being  examined. 

With  the  area  inside  the  horizontal  triangle  formed  by  the  three  bal¬ 
loons  computed  for  each  minute  using  (2),  equation  (1)  was  reviritten  as 


- _ i _ _  iki±L  ~  n't 

(At +  Vi)  6c> 


(b) 


where  t  is  again  the  time  of  reading  of  the  theodolites  and  with  t  ■  0  at 

release  time.  The  value  of  the  in  sec-1,  from  (h)  then  represents 

H  H 

the  mean  divergence  between  levels  determined  by  t  and  (t+1). 

The  vertical  component  of  the  relative  vorticity  is  given  by 


We  can  define  a  vector  tyV  (fig.  2)  which  is  equal  in  magnitude  to  the  ve¬ 
locity  vector  but  is  directed  90*  to  the  right  of  [2j.  The  two- 
dinensiona-d  divergence  of  this  vector  is  then, 

w-fe-jy-s 

The  coord i nates  at  the  tine 
t  of  the  fictitious  position  of  the 
balloon  n  (where  n  is,  as  before,  1, 

2,  or  3>  depending  on  the  release  site 
under  consideration)  which  moves  sub¬ 
ject  to  the  "wind"  from  observa¬ 
tion  site  n  are  given  by 


(6) 
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Fig.  2.  Definition  of  vector 
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v  •  «  :  .  v  *  m  V 

no  no  *  J  no  ‘'no 

UXl  ( x'nk  "  KA-l  +  ^Jt-1 

t  (7) 

=  (y‘n}t-l  -  ^XH-1 

for  t  ■  1,  2,  3>  ...  Note  that  at  release  time  the  initial  coordinates 
of  and  at  sites  1,  2,  and  3  are  identical.  Therefore,  in  equa¬ 
tions  7,  the  subscript  t  cannot  be  equal  to  zero,  so  that  in  these  equa¬ 
tions  t  represents  the  tine  of  the  reading,  in  minutes,  after  release. 

*  ..  and  refer  to  the  differences  from  one  reading  to  the  next  of 
the  coordinates  of  vector as  determined  by  equations  3. 

The  vorticity  of  is  exactly  the  divergence  ofW^{  ana  so  is 
given  by  relations  similar  to  2  and  1, 

x'i  y'i  1 

*'t  -  1/2  2  y«2  1  (&) 

x* 3  y*3  1 


where  all  quantities  have  a  like  meaning,  except  that  the  primed  quanti¬ 
ties  refer  to  vector  the  rotated  vector. 

3.  Sample  calculation 

The  method  used  to  determine  the  divergence  is  perhaps  best  shown 
by  an  actual  example.  The  three  minute  ascension  of  8  June  19?!j  (l?OC  2) 
was  chosen  for  this  purpose. 

Table  1  shows  a  sample  work  sheet  for  the  computation  of  the  x  and  y 
coordinates  for  the  balloons  released  at  sites  1,  2,  and  3.  Data  in  paren¬ 
theses  in  this  table  are  the  raw  data  as  recorded  by  the  observing  teams. 
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Table  1.  Sample  worksheet  for  the  computation  of  the  x  anti  y  coordinates 
of  the  three  balloons  at  the  end  of  each  minute  interval.  (Release: 
1000  LIST }  8  June  19*6.) 


zzzzzzzzzzzzzzzzzzzzzzzzazzzz, 

Site  1 


ilinute 

(ex) 
cot  e^ 

r  R1 

[feet j 

(0(l) 

sin  oC^ 

r  X1  , 
[feet] 

(oti; 

cos  V 

yi 

[feet] 

1 

(32.1(8 ) 
1.5707 

1093 

(297.56) 

-0. 06647 

5852 

(297.56) 

0.46278 

1763 

2 

(30.31) 

1.3605 

1846 

(304.00) 

-0.82904 

6413 

(304.00) 

0.55919 

1257 

3 

(34.82) 

1.6379 

2859 

(3o8.55) 

-0.78208 

7115 

(308.55) 

0.62320 

507 

Site  2 

Minute 

(e2) 

cot  _0 

c. 

f  R2 
(feet) 

(<*2) 
sin  ^2 

x2 

[feet j 

(*2) 
cos  oi2 

y2 

[feet] 

1 

(52.16) 

0.77661 

541 

(316.31) 

-C. 69067 

374 

(316.31) 

0.72317 

5557 

2 

(51.83) 

0.78596 

1067 

(312.38) 

-0.73865 

788 

(312.38) 

0.67409 

5229 

3 

(49.29) 

0.86064 

1711 

(313.80) 

-0.72176 

1235 

(313.80) 

0.69214 

4764 

Site  3 

Minute 

(23) 
cot  e^ 

R, 

[feet] 

(<X3) 

sin  ^3 

x3 

[feet] 

(*3) 

COS  "*(-3 

r  y3  - 
[feetj 

1 

(53.32) 

0.74492 

516 

(302.35) 

-0.84480 

2292 

(302.35) 

0.53509 

-277 

2 

(54.91) 

0.70238 

953 

(298.75) 

-0.87673 

2690 

(298.75) 

0.68099 

-658 

3 

(49.35) 

0.85862 

1707 

(309.62) 

-0.77033 

3169 

(309.62) 

0.63765 

-1088 
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Table  2,  Divergence  computation  worksheet-  (Release:  1000  EST, 
8  June  1951* - ) 


Time  (minutes) 

Coordinates  0  12  3 

(feet) 


k 

1 

ft 

•« 

■i 

•  *ik 

:*«■*(»«*'*  wta 

x,  1*663 

(4xx) 

(969) 

5852 

(561) 

61*13 

(706) 

7119 

Xp  0 

(^x2) 

07k) 

371* 

(1*11*) 

786 

(1*1.7) 

1235 

X.  18  51* 

"  (4x3) 

(1*36) 

2292 

(398) 

2690 

(1*79) 

3169 

y,  2289 

1  (£yx) 

(-506) 

1763 

(-526) 

1257 

(-750) 

507 

y..  59i*8 

2  (<$y2) 

(-291) 

5557 

(-328) 

5229 

(-1*65) 

1*761* 

y->  o 

"  (^y,) 

(-277) 

-277 

(-181) 

-1*58 

(-630) 

-1CB8 

Area  11,130,11*9  12 

(sq.  ft.) 

,360,0fc0  12 

,217,316 

] 

13, 100,065 

Divergence 
(xl0'5  sec-1 ) 

17U.53 

-19.36 

116.22 

The  coordinates 

of  the  balloons 

released 

at  sites 

1,  2, 

and  3  are 

recorded  again  in  table  2.  The  quantities  in  parentheses  in  this  table 
are  the  differences  between  the  coordinates  from  one-minute  reading  to 
the  next  as  inaicateu  in  the  left  hand  column.  The  lower  part  of  the 
table  gives  the  area  at  each  minute,  from  equation  2,  and  the  mean  diver¬ 
gence  between  each  minute  reading  from  equation  I*. 

Table  3,  which  is  similar  to  table  2,  shows  the  coordinates  of  the 
three  rotatea  vectors  of  sites  1,  2,  and  3  as  computed  from  the  values  in 
parentheses  in  table  2  and  equations  7.  The  areas  of  the  fictitious 


7 


-<’l 

1*883 

1*377 

3851 

3101 

xs 

Cm 

0 

-391 

-719 

-1181* 

Xb 

1351* 

1577 

1396 

76c 

2289 

1  j20 

759 

...| 

591*8 

557U 

5160 

1*713 

,1 

''  3 

0 

-);33 

-836 

-1315 

Ar*af 
(sq,  ft.) 

11,130,11*9 

10,11)6,670 

9,0/4  6,803 

8,371,1*90 

Vorticity 

(::10“5  scc-1)  -15U.08  -191.01  -129.23 


triangles  formed  by  the  rotation  of  the  three  horizontal  range  vectors 
at  each  minute  reading,  from  equation  8,  are  shown  in  the  bottom  portion 
of  the  table,  along  with  the  vorticity  values  computed  from  equation  9. 

It,  The  observations 

As  in  the  study  of  the  ascentional  rates  [lj,  the  three  balloons 
for  each  run  were  filled  in  the  Weather  Bureau  helium  shelter.  In  this 
case,  however,  it  was  not  possible  to  check  for  gas  leakage  as  carefully 
as  for  the  ascensional  rate  check  releases,  but  again  the  balloons  were 
filled  to  the  correct  free  lift  as  carefully  as  possible  with  the  standard 
equipment.  Simultaneous  release  of  the  three  balloons  was  accomplished 
as  before  by  stopwatches  which  were  preset  before  taking  the  balloons  and 
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equipment  to  the  release  sites.  Again,  all  balloon  ascensions  were  made 
by  two-nuan  observation  teams  reading  both  the  azimuth  and  elevation  angles 
to  0.01  degrees. 

Altogether  fourteen  successful  triple-balloon  runs  were  accomplished, 
seven  day  runs  and  seven  night  runs.  The  day  runs  averaged  about  28,000  ft 
(excluding  the  three-minute  run  of  June  8)  and  the  night  runs  averaged  about 
16,000  ft. 

The  data  obtained  were  set  up  on  punched  cards  for  machine  computa¬ 
tion  of  the  divergence  and  vorticity.  Any  of  the  machine  computed  results 
which  appeared  to  be  out  of  line  were  checked  by  hand  and  corrected  when- 
even  necessary.  It  is  felt  that  the  resulting  values  of  divergence  and 
vorticity  presented  herein  are  free  of  computational  error  to  a  high  de¬ 
gree.  There  is  one  possible  exception,  and  this  particular  run  (on  June  6 
at  1500  Z)  will  be  discussed  in  some  detail  in  section  6.  It  has  been 
omitted  from  the  divergence  and  vorticity  mean  values  discussed  below. 

5.  Analysis  of  the  data 

All  of  the  divergence  and  vorticity  values  obtained  are  based  on 
one^ninute  changes  in  area  and  orientation  of  a  small  triangle  (about 
one  mile  on  a  side)  as  this  triangle  moved  through  a  vertical  distance 
of  about  620  ft.  The  detail  of  the  small  scale  wind  field  was  such 
that  some  smoothing  of  the  data  was  desirable  in  order  to  see  the  main 
features  revealed  by  the  study. 

a.  The  divergence 

The  divergence  computations  for  each  run  were  smoothed  by  taking 
five-minute  means  in  the  vertical  which  had  a  two-minute  overlap.  That 
is,  a  mean  value  was  obtained  at  three-minute  intervals  (3,  6,  9,  etc) 
where  the  sixth-minute  value  was  an  average  of  the  fourth  through  eighth 


9 


minutes,  etc.  Then  the  averages  at  each  three-minute  interval  for  all  day 
runs,  all  night  runs  and  all  runs  (day  and  night)  were  computed. 

Mean  curves  of  day  runs,  night  runs  and  all  runs  versus  height  are 
shown  in  fig.  3.  It  is  interesting  to  note  that  the  0300  Z  and  1500  Z 
curves  are  quite  similar  to  about  15,000  ft  which  is  about  the  extent  of 
the  0300  Z  data.  The  last  plotted  point  on  this  curve  (twenty -seventh 
minute)  is  a  mean  of  only  three  runs  while  the  point  just  below  (twenty- 
fourth  minute)  is  an  average  of  five  runs.  Also,  the  last  four  points  on 
the  1500  Z  curve  are  means  of  no  more  than  three  or  two  runs  while  the 
fourty-second  minute  point  is  a  mean  of  five  runs.  Both  curves  are  dashed 
in  the  region  where  a  sample  smaller  than  five  runs  is  involved.  In  the 
0300  Z  curve  there  is  a  decrease  in  divergence  with  height  reaching  a 
minimum  at  about  6000  ft,  then  increasing  to  a  maximum  just  below  12,000 
ft,  showing  a  secondary  minimum  at  about  13,500  ft  and  then  increasing 
again  (the  top  point  must  be  considered  doubtful).  The  1500  Z  curve  shows 
the  two  minima  to  be  about  1500  ft  hi jher  than  they  were  found  in  the  0300 
Z  curve  but  the  maximum  appears  again  just  below  12,C0U  ft. 

At  higher  levels  the  1500  Z  curve  increases  to  a  maximum  at  about 
22,000  ft,  decreases  to  another  minimum  at  about  26,000  ft  (again  con¬ 
sidering  as  doubtful  the  points  based  on  three  runs  or  less)  and  appears 
to  increase  again  at  30,000  ft  and  above. 

Returning  to  the  lowest  15,000  ft  of  the  two  curves  and  the  overall 
mean  curve,  there  are  some  noteworthy  similarities  and  differences  when 
compared  with  a  similar  study  of  larger  scale  motion  [3j.  In  [3]  it  was 
found  that  the  mean  divergence  over  an  area  of  about  2  x  10^  niles^  showed 
some  interesting  diurnal  variations,  namely,  the  divergence  at  all  levels 
(below  20,0C0  ft)  was  greater  at  0300  Z  than  1500  Z  and  the  levels  of 


(U  P  spsL|4)  m&!3H 
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least  divergence  were  lower  at  1500  Z  than  at  0300  Z.  The  mean  0300  Z 
and  1500  2  curves  in  this  study  seem  to  be  contradictory  to  the  results 
in  [3J  in  that  0300  Z  has  smaller  divergence  values  at  all  levels  than 
1500  2  and  the  levels  of  least  divergence  are  lower  at  0300  Z  than  at 
1500  2%  However,  the  area  in  the  present  study  (**£  mile2)  is  only  a 
point  when  compared  with  the  area  in  [3].  There  is  no  reason  to  suspect 
that  every  point  in  the  area  considered  in  13]  would  give  the  same  re¬ 
sult  as  the  total  area.  In  fact,  it  would  be  surprising  if  such  were  the 
case. 

Tne  similarities  with  [3J  are  that  the  level  of  least  divergence 
occurs  at  6000  ft  in  one  case  and  8000  ft  in  the  other  and  that  a  secon¬ 
dary  level  of  minimum  divergence  occurs  at  about  13,000  ft  in  both  cases 
with  an  apparent  increase  above  this  level.  It  may  be  quite  significant 
that  in  two  cases  which  deal  with  different  seasons  (spring  and  summer), 
different  regions  (both  in  size  and  location)  end  vastly  different  scales 
of  motion,  that  the  level  where  the  divergence  is  noticeably  least  is  in 
the  6000  -  30C0  ft  region.  This  is  rather  low  for  a  level  of  non-divergence 
as  it  is  generally  thought  of.  Also  it  was  shown  in  [bj  that  the  level 
of  least  divergence  found  in  [3 J  wa s  the  level  where  absolute  vorticity 
was  best  advected  with  observed  winds. 

Fig.  U  may  shed  some  light  on  the  difference  in  magnitude  in  the 
mean  divergence  at  0300  2  and  1500  Z.  This  is  a  curve  connecting  points 
of  mean  divergence  for  each  run  averaged  over  all  heights.  The  mean  diver¬ 
gence  was  apparently  larger  in  early  June  than  in  late  June  and  early 
July,  regardless  of  time  of  day.  The  fact  that  the  early  June  runs  were 
nearly  all  at  1500  Z  may  be  the  sole  reason  that  the  mean  1.50O  Z  curve 
shows  larger  divergences  level  for  level  than  the  mean  0300  Z  curve. 


0300Z 
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Table  ij  shows  the  mean  values  of  divergence  for  day  runs,  night  runs 
and  all  runs.  The  sign  of  the  divergence  was  disregarded  in  obtaining  the 
means.  The  data  in  table  l^are  confined  to  the  lowest  17, OCX)  ft  of  the 
atmosphere.  Five  of  the  seven  0300  Z  runs  went  as  high  as  twenty-seven 
minutes  17,000  ft)  but  only  three  went  higher,  two  of  these  three 
terminating  at  the  twenty-ninth  minute.  Host  of  the  1500  Z  runs  went  con¬ 
siderably  higher  but  since  the  0300  Z  and  1500  Z  times  are  averaged  to¬ 
gether  to  get  an  overall  value^a  common  thickness  must  be  used.  Therefore, 
only  the  first  twenty-seven  minutes  (**  17,000  ft)  are  used  in  determining 
the  table  k  figures. 


Table  h 

1500  Z 
(day  runs) 

0300  Z 
(night  runs) 

All  runs 

Number  of 
computations 

13b 

160 

298 

Mean  divergence 
:c  10"5  sec-1 

ob.  9 

h2.o 

5*1.8 

The  overall 

time  and  height  mean  divergence  value  of  5*1.8 

;c  10“5  sec” 

is  a  little  less  than  one  and  one-half  orders  of  magnitude  greater  than 
the  overall  mean  value  of  1.5  x  10“^  sec"*  found  for  the  much  larger  scale 
motion  in  [3]. 

The  above  divergence  discussion  has  dealt  entirely  with  smoothed  curves 
and  mean  values.  As  mentioned  earlier  the  individual  runs,  especially  when 
the  sign  is  considered,  show  much  detail  of  small-scale  atmospheric  motions. 
An  individual  run  showing  unsmoothed  curves  of  divergence  and  vorticity 
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versus  height  is  presented  in  fig.  6  (section  6).  It  is  felt  that  the  many 
oscillations  from  plus  to  minus  and  vice  versa  are  factual  (containing,  of 
course,  a  percentage  of  error  which  will  be  discussed  below).  However,  we 
cannot  compare  these  small  scale  oscillations  directly  with  other  studies 
(such  as  [3])  or  with  our  synoptic  maps  due  to  the  vast  difference  in  scale 
of  the  observational  networks  involved.  Therefore,  the  small  oscillations 
were  deliberately  smoothed  out  so  that  such  comparisons  could  be  made.  Some 
details  of  fig.  6  will  be  discussed  in  section  6. 
b.  The  vorticity 

The  same  smoothing  technique  was  applied  to  the  vorticity  as  was 
mentioned  in  the  divergence  discussion  above.  Fig.  5  shows  the  mean  curves 
for  0300  Z,  1500  Z  and  all  runs  averaged  together.  As  in  the  divergence, 
the  0300  Z  vorticity  is  less  at  all  levels  than  the  1500  Z  vorticity.  Also 
the  three  curves  show  quite  a  bit  of  similarity  to  the  divergence  curves, 
having  maxima  and  minima  at  about  the  same  levels  up  to  17,000  ft. 

Table  5  shows  mean  vorticity  values  for  day  runs,  night  runs  and 
all  runs.  The  overall  time  and  height  mean  of  vorticity  *  7I4.I4  x  10~5  sec”* 
is  about  I.I4  times  as  great  as  the  similar  value  for  the  divergence  in 


Table  5 


1500  Z 

day  runs 


0300  Z 
night  runs 


All  runs 


Number  of 


computations 

138 

160 

298 

Mean  vorticity 
x  10“5  sec"* 

108.9 

UU.6 

7h.h 

is 


Cw  JO  spsui)  |q6i8H 


sec 
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table  U.  Again  only  the  first  twenty-seven  minutes  of  data  are  included  in 
the  mean  values  in  table  t>. 

The  mean  vorticity  values  for  each  run  showed  the  sane  general 
trend  evident  in  the  divergence  in  fig.  lu  That  is,  the  larger  mean  vor- 
ticities  occurred  early  in  June  with  smaller  values  in  late  June  and  early 
July.  Again  the  fact  that  the  early  June  runs  were  nearly  all  1^00  Z  runs 
could  be  the  principal  reason  for  the  apparent  diurnal  difference  seen  in 
fig.  Also  the  vorticity  shows  somewhat  greater  percentual  differences 
between  the  0300  Z  and  1500  Z  mean  curves  (fig.  5)  than  does  the  diver¬ 
gence  (fig.  3). 

Although  the  mean  vorticity  and  divergence  curves  in  figs.  5 
and  3  appear  to  indicate  a  strong  correspondence  in  maximum  and  minimum 
values  as  far  as  elevation  is  concerned,  the  individual  runs  did  not  show, 
minute  for  minute,  any  recognizable  correspondence  of  this  nature  even  when 
absolute  values  were  compared.  So  that  the  possible  implication  that  the 
correspondence  of  mean  values  night  be  inherent  in  the  computing  method 
is  not  valid  in  this  case. 

6.  The  triple-balloon  run  of  June  6,  1 95h  at  1500  Z 

The  June  6th  run  showed  some  unusually  large  divergence  values  (as 
large  as  +  2  x  10"^  sec"*)  and  was  therefore  excluded  from  the  averages 
shown  in  table  h  and  fig.  3.  a  more  detailed  investigation  of  this  run 
serves  to  answer  some  questions  concerning  it  but  fails  to  prove  conclu¬ 
sively  whether  these  large  divergence  values  are,  or  are  not,  real. 

The  very  large  divergence  values  occurred  in  the  layer  from  the 
twenty-second  to  thirty-second  minutes  of  ascent  (approximately  13,000 
to  19,000  ft).  The  curves  of  divergence  and  vorticity  versus  height  are 
shown  in  fig.  6.  It  should  be  mentioned  that  the  assumed  rate  of  ascent. 


o  m  o  m  o  m  O 

ro  ro  c\l  c\j  —  — 
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rather  tiian  the  day  rate  determined  in  [1]  was  used  in  working  up  the  three 
runs  so  that  direct  comparison  with  winds  at  other  stations  would  be  possible. 
The  D  and  values  however  are  based  on  a  rate  [1]  which  placed  the  balloons 
about  600  ft  (or  one  minute)  higher  in  the  twenty-second  to  thirty-second 
minute  region.  It  will  be  noted  that  the  vorticity  values  in  the  13,000 
to  19,000  ft  layer  are  of  the  same  magnitude  as  at  all  other  levels.  This 
is  because  the  triangle  formed  by  the  balloons  was  experiencing  very  little 
rotation  but  its  area  was  changing  very  rapidly. 

When  this  run  was  taken,  the  balloon  at  site  2  was  erroneously  released 
two  minutes  before  the  prearranged  release  time.  Immediately  recognizing 
the  mistake,  the  observing  teem  decided  to  follow  the  site  1  balloon  to 
get  a  rate  of  ascent  check  (see  [lj).  Later,  realizing  that  the  tvo-minute 
difference  could  be  adjusted,  they  switches  to  the  site  2  balloon  at  the 
seventh  minute  and  recorded  it  as  the  ninth  minute  since  this  balloon  had 
been  released  two  minutes  early.  Although  the  height  of  the  site  2  balloon 
run  was  corrected  for  the  time  difference,  it  might  be  suspected  that  this 
difference  in  time  was  instrumental  in  causing  the  large  divergences  which 
were  computed,  examination  of  the  vind  direction  versus  height  reveals, 
however,  that  sites  2  and  3  are  almost  identical,  particularly  from  the 
seventeenth  minute  (10,000  ft)  upward  (see  fig.  10).  It  was  the  site  1 
ascent  which  differed  from  the  other  two  ascents  in  the  layer  of  large  diver¬ 
gence  values.  A  pronounced  change  in  wind  direction  (from  2NE  to  W  through 
N)  took  place  in  the  site  2  and  3  data  between  the  twenty-second  and  twenty- 
fourth  minutes,  but  this  same  change  did  not  take  place  in  the  site  1  data 
until  the  twenty-fourth  to  twenty-sixth  minutes.  The  wind  speeds  versus 
height  bear  out  the  direction  curves,  with  the  speed  minimum  at  each  site 
occurring  at  the  same  height  as  the  sharp  change  in  wind  direction. 
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The  effect  of  the  rapid  wind  uirection  change  acting  two  Minutes  later 
at  site  1  than  at  sites  2  and  3  is  illustrated  in  fig.  7.  The  position  and 
orientation  of  the  triangle  is  shown  at  tine  zero  (dashed  triangle)  and  at 

i 

various  minute  readings  during  the  ascent.  Throughout  the  layer  shown 
(twenty-one  to  thirty-three  minutes  of  ascent),  the  side  of  the  triangle 
connecting  sites  2  and  3  remained  almost  exactly  parallel  with  its  original 
orientation.  The  triangle  of  balloons  was  moving  toward  the  US'J  until  the 
twenty-third  minute  at  which  time  the  site  2  and  3  balloons  started  moving 
toward  the  ESd.  The  site  1  balloon,  however,  continued  to  move  WSJ  until 
the  twenty-fifth  minute  and  this  caused  the  area  of  the  triangle  to  decrease 
rapidly,  actually  going  to  zero  just  before  the  twenty-sixth  minute  reading. 
The  area  shows  an  increase  again  to  the  twenty-ninth  minute,  a  decrease  from 
the  twenty-ninth  to  thirty-first  minutes  and  then  another  increase.  The 
divergence  curve  in  fig.  6  has  been  broken  to  indicate  that  theoretically 
it  goes  to  and  returns  from  infinity  as  the  triangular  area  goes  to  zero. 

It  can  be  seen  that  the  amount  of  turning  of  the  triangle  was  small,  there¬ 
fore  the  vorticity  values  in  the  layer  were  quite  small,  averaging  only 
9U  x  1CT^  sec”*. 

There  are  at  least  two  possible  explanations  of  the  observed  results; 
a)  the  site  1  balloon  could  for  some  reason  be  lower  than  the  other  two  by 
about  1200  ft,  and  therefore  not  experiencing  the  direction  change  at  the 
same  time,  or  b)  the  site  1  balloon  could  be  actually  experiencing  the  wind 
direction  change  at  a  higher  level  than  the  other  two  balloons.  In  connec¬ 
tion  with  a),  the  idea  immediately  occurs  that  the  site  1  balloon  had  a 
pinhole  and  was  leaking  gas.  However,  all  balloons  were  checked  for  pin¬ 
holes  when  inflated.  Also  after  the  balloons  were  inflated  a  period  of 
at  least  fifteen  minutes  passed  before  they  were  released,  due  to  the  fact 
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that  all  balloons  were  filled  at  site  1  and  then  two  ol'  them  were  trans¬ 
mitted  to  sites  2  ana  3*  bven  if  a  pinhole  had  gone  unnoticed,  this  time 
lag  would  Uav:  caused  sufficient  gas  loss  before  release  so  that  the  leaky 
balloon  would  rise  more  slowly  throughout  the  entire  run  with  the  difference 
in  rates  becoming  or:  no  tic.- bib  at  higher  levels.  In  fact,  in  one  of  the 
runs  it  was  felt  that  exactly  this  had  happened.  The  direction  curves  versus 
height  of  two  sites  were  nearly  identical  with  the  third  showing  the  same 
features  but  at  different  levels.  In  this  case  the  lag  was  one  minute  in 
the  surface  to  ninth  minute  layer,  two  rinutes  in  the  tenth  to  twenty-fourth 
minute  layer  end  three  minutes  in  the  twenty-fifth  to  thirty-ninth  minute 
layer.  However,  in  the  June  l,  1500  Z  run  there  is  a  remarkable  correspon¬ 
dence  betv/cen  the  cirection  curves  of  the  three  balloons  up  to  about  the 
twenty-second  minute  13,uOO  ft),  then  a  two  minute  uifference  appears  in 
the  direction  changes  occurring  in  the  twenty-second  to  thirty-second  minute 
layer  (^-  13,000  to  15,000  ft).  Above  15,000  ft  the  direction  changes  with 
height  are  so  slight  that  no  correspondence  or  contrast  can  be  noted.  This 
would  make  the  possibility  of  a  preexisting  pinhole  seen  unlikely.  Also,  if 
a  leak  had  in  some  way  developed  at  alxrnt  the  twenty-first  or  twenty-second 
minute  and  was  large  enough  to  cause  a  two-minute  lag  almost  immediately, 
then  one  would  expect  the  lag  to  increase  rapidly  as  time  went  on.  But  the 
lag  did  not  increase  noticeably  as  shown  by  the  plotted  directions  in  the  layer 
between  the  twenty-sixth  and  thirty-second  minutes  (see  fig.  10). 

It  hardly  seems  plausible  that  the  site  1  balloon  is  some  1200  ft  lower 
than  the  other  two  in  the  layer  under  discussion  due  to  differential  vertical 
motions  of  a  small  scale.  There  is  pretty  yeneral  agreement  among  researchers 
that  +  10  - —  is  a  large  vertical  velocity  value,  and  this  is  only  3  per  cent 

m m  S  00 

of  the  rate  of  ascent,  even  if  larger  vertical  velocities  are  assumed  for 
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smaller  scale  motions  such  as  may  be  measured  here,  it  is  still  difficult 
to  visualize  one  balloon  falling  1200  ft  behind  the  other  two  in  a  period 
of  about  one  or  two  minutes  due  to  differential  vertical  velocities.  The 
only  clouds  present  during  this  ascent  were  scattered  cirrus  at  about 

30,0-00  ft. 

Now,  there  is  also  the  possibility  (b)  that  the  pronounced  change  in 
wind  direction  actually  occurred  at  a  higher  level  at  site  1  than  at  sites 
2  and  3.  It  is  largely  a  product  of  our  thinking  in  terms  of  the  existing 
upper  wind  network  that  we  would  accept  as  normal  a  130°  change  of  wind 
direction  in  the  vertical  occurring  say  at  lU,000  ft  at  one  station  and  at 
16,000  ft  at  a  station  100  to  200  miles  distant,  but  would  be  quite  reluc¬ 
tant  to  accept  the  idea  that  this  change  might  occur  at  lii,000  ft  at  one 
station  and  at  15,000  ft  at  another  if  the  stations  were  only  one  mile  apart. 
There  is  no  reason  why  this  night  not  happen  but  we  are  not  accustomed  to 
thinking  that  it  does  because  we  do  not  normally  see  it  taking  place. 

A  look  at  the  upper  wind  pattern  over  the  extreme  southeastern  United 
States  at  this  time  revealed  that  an  anticyclone  existed  at  all  levels  up 
to  18,000  ft  in  southern  Alabama.  Detailed  isogon-isovel  analyses  of  the 
levels  12,000,  lit, 000,  16,000  and  18,000  ft  indicated  that  the  anticyclone 
was  displaced  slightly  westward  at  each  succeeding  level  (see  fig.  8).  The 
companion  hyperbolic  point  was  found  south  of  Jacksonville  at  12,000  ft  and 
was  displaced  westward  to  Pensacola  at  18,000  ft.  Fig.  8  shows  the  positions 
of  the  anticyclone  (A)  and  the  hyperbolic  point  (H)  at  the  levels  mentioned 
above.  At  20,000  ft  the  point  singularities  could  not  be  found,  suggesting 
that  the  analysis  at  that  level  should  show  a  ridge  open  to  the  south.  This 
is  reasonable  when  one  considers  the  trend  of  (A)  and  (H)  with  height  shown 
in  fig.  8. 
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Streamlines  drawn  from  ruled  isogons  and  from  plotted  wind  vectors 
both  emphasized  the  elongated  or  elliptical  shape  of  this  anticyclone,  the 
major  axis  running  WNW-ESE.  A  freehand  streamline  analysis  at  lli,000  ft 
is  shown  in  fig.  9.  Here  the  elongated  shape  of  the  anticyclone  has  been 
emphasized  more  than  in  the  customary  analysis  which  is  always  smoothed 
to  a  scale  consistent  with  the  distance  between  wind  reports.  That  is, 
instead  of  assuming  that  the  wind  direction  changes  in  a  regular  fashion 
from  the  WMVf  reports  at  Maxwell  AFB,  Birmingham  and  Atlanta  to  NE  at 
Tallahassee  (some  200  to  250  miles  away),  it  has  been  assumed  that  the 
WNW  flow  continues  to  within  a  short  distance  of  Tallahassee.  If  this 
short  distance  is  taken  to  be  about  one-half  mile,  then  the  directions 
observed  at  lb, 000  ft  at  the  theodolite  sites  2  and  3  are  consistent  with 
the  larger  scale  pattern.  The  direction  change  had  nearly  all  taken  place 
at  sites  2  and  3  by  the  time  the  balloons  had  reached  lb, 000  ft.  The  site 
1  balloon  direction  change  was  conpleted  just  above  15,000  ft  so  that  it 
still  shoved  a  northeasterly  wind  at  lb, 000  ft.  It  can  be  seen  in  fig.  8 
that  the  hyperbolic  point  was  about  at  the  longitude  of  Tallahassee  at 
15,000  ft.  As  soon  as  this  point  assumes  a  position  west  of  the  station 
(TLH)  the  station  wind  must  begin  to  have  a  westerly  component.  Although 
the  westward  motion  of  the  hyperbolic  point  can  explain  the  shift  in  wind 
direction  from  ENE  to  WNW  at  the  station  (station  here  refers  to  site  1, 
the  official  Weather  Bureau  wind),  it  does  not,  or  course,  explain  the  direc¬ 
tion  shift  at  sites  2  and  3  which  occurred  mostly  below  lit, 000  ft.  How¬ 
ever,  if  one  considers  that  a  sharp  line  of  shear  may  exist  between  north- 
westerlies  and  northeasterlies  and  that  this  line  lies  between  site  1  and 
sites  2  and  3  at  lb, 000  ft,  then  it  is  simple  to  show  (fig.  10)  how  the 
observed  distribution  of  three  local  northeast  winds  at  12,000  ft,  two 
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northwest  winds  ami  one  northeast  xjind  at  1)4,000,  and  tliree  local  northwest 
winds  at  16,000  ft  and  higher  could  occur.  The  heavy  line  in  fig.  10  shows 
the  separation  of  easterly  from  westerly  winds.  We  are  accustomed  to  think¬ 
ing  of  this  heavy  line  as  nearly  horizontal  over  such  small  linear  distances 
as  3029  ft,  which  is  the  east-west  distance  between  sites  3  and  1.  In  fig. 
10,  where  there  is  a  one  to  one  ratio  between  ordinate  and  abscissa,  a 
nuxiMum  tilt  of  21*  from  the  horizontal  would  be  required  for  the  heavy 
line  and  this  maximum  would  occur  between  sites  3  and  1. 

This  argument  cannot  be  definitely  proved  nor  disproved  from  the  data 
available.  It  is  an  entirely  possible  e.qjlanation  which  allows  the  data 
from  the  three  sites  to  be  consistent  locally  and  consistent  with  the  larger 
scale  pattern.  Those  working  with  wind  fields  on  a  small  scale  are  contin¬ 
ually  fin  ing  evidence  that  such  things  way  be  more  co.-aion  than  generally 
thought.  ..cLeen  and  kedos  [£]  recently  noted  a  very  sharp  shear-line  at 
36, too  ft.  It  was  not<a>  from  aircraft  observed  winds  anti  could  not  be  found 
with  the  coverage  of  conventional  winds.  This  shear-line  (located  about 
110  miles  cast  of  Tallahassee  in  April  1 9$h)  showed  a  18G*  wind  direction 
change  in  a  distance  of  about  2b  riles  (the  distance  between  observations) 
with  s,  cebs  of  ;fxmt  ?0  knots  on  each  sice  of  the  line  of  direction  change. 
This  she;  r-line  night  have  shown  up  as  even  sharper  if  observations  had  been 
available  ;t  smaller  'distance  intervals.  The  ^cLean-uad.os  shear  -line  w;s 
fount  on  c;ch  of  two  successive  passes  by  the  aircraft  through  the  region 
at  times  about  two  hours  apart,  except  for  minor  details  this  shear-line 
fit  into  a  conventional  contour  analysis  at  the  sane  level  and  time}  how¬ 
ever  much  additional  information  \^as  given  by  the  aircraft  observations. 

7.  discussion  of  errors 

An  attempt  has  been  male  to  determine  as  nearly  as  possible  the  per- 
centunl  error  included  in  the  computed  divergence  and  vorticity  values. 
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Tlie  three  largest  sources  of  error  are: 

1)  Whatever  difference  exists  in  rates  of  ascent  of  three  balloons 
released  simultaneously  in  a  small  region,  2)  the  difference  between  the 
rate  of  ascent  on  any  given  day  or  night  and  the  mean  ascent  rates  used 
for  the  computations  and  3)  the  errors  in  the  azimuth  angle  readings  - 
particularly  when  the  flight  of  the  balloon  causes  rapid  time  changes  in 
azimuth. 

Error  1:  The  error  caused  by  differences  in  ascent  rates  of  three  bal¬ 
loons  released  simultaneously  is  not  available  from  the  data 
gathered.  However  such  differences  as  do  exist  would  be  due 
to  very  snail  scale  atmospheric  vertical  motions,  since  the 
balloons  were  released  in  a  small  region,  and  should  be  quite 
random . 

A  spot  check  was  made  of  divergence  values  of  (10,  50,  100, 

150,  200,  300)  x  1C~£  sec”1  by  introducing  an  assumed  error 
in  heights.  At  one  minute  the  heights  originally  used  were 
changed  by  (+10,  0,  -10)  ft  and  heights  at  the  next  minute  were 
left  unchanged.  This  amounts  to  including  a  20  ft  error  be¬ 
tween  stations.  The  divergence  values  were  recomputed  and  the 
percentual  errors  are  shown  in  table  6.  Two  such  spot  checks 
were  made  for  each  divergence  value  listed  above  and  the  check 
points  included  nearly  all  of  the  triple-balloon  runs, 

Error  2:  It  was  found  that  the  individual  runs  in  Part  1  differed  from 

the  mean  rates  of  ascent  by  18  ft/min  on  the  average,  or  a  per- 
centual  difference  of  <  3  per  cent  of  the  mean  ascent  rate.  It 
is  difficult  to  tell  how  much  the  computations  are  in  error  due 
to  this  effect.  It  is  probable  that,  on  a  given  day  or  night, 
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all  three  balloons  will  rise  faster  or  slower  or  the  same  as 
the  mean  rate  used  in  the  computations.  It  has  been  assumed 
that  the  one  minute  thicknesses  are  18  ft  in  error  due  to  this 
effect  and  the  same  percentual  figures  (see  table  6)  have  been 
used  that  were  applied  in  the  case  of  error  1  where  a  20  ft 
error  was  assumed. 

error  3:  At  the  same  check  points  used  in  determining  error  1,  azimuth 

angle  reading  errors  of  (+.02,  0,  -.02)  degrees  were  introduced 
into  the  original  angles  at  the  first  minute  and  the  readings 
at  the  next  minute  were  left  unchanged.  The  divergence  values 
were  recomputed  and  the  results  are  shown  in  table  6. 
irror  Us  This  error  has  been  added  to  the  table  to  take  care  of  the  various 
very  small  errors  such  ass 

a.  A  slight  error  which  is  introduced  due  to  human  inability 
to  exactly  orient  the  three  theodolites,  thereby  causing  a  very 
slight  random  tilt  in  the  original  triangle. 

b.  The  ground  was  not  exactly  flat  and  this  introduces  a 
slight  systematic  tilt  in  the  original  triangle.  Both  a  and  b 
are  extremely  small  and  are  carried  through  the  entire  run.  b 
will  be  the  seeae  for  all  runs.  These  two  errors  are  considered 
negligible. 

c.  The  computing  method  is  such  that  very  snail  angles  in 
the  geometry  of  the  setup  can  cause  large  errors.  It  is  felt  that 
this  error  lias  been  reduced  to  the  negligible  class  by  eliminating 
all  computations  which  contained  trigonometric  function  values 
smaller  than  or  equal  to  0.10000, 

d.  Readings  of  angles  may  not  have  been  made  at  exactly  the 


same  time  at  the  three  installations.  In  all  cases  the  readings 


30 


were  either  on  time  or  late,  never  early.  In  all  cases  of  late 
readings  an  attempt  was  made  to  note  the  lateness  so  that  the 
values  could  be  suitably  adjusted  later  when  the  data  were  worked 
up.  It  is  felt  that  this  procedure  reduced  this  error  to  the 
negligible  class. 

e.  Only  brief  mention  will  be  given  to  errors  in  reading 
elevation  angles  because;  1)  this  angle  changes  slowly  and  can 
be  read  more  accurately  than  the  azimuth,  and  2)  any  error  which 
might  exist  here  is  not  independent  of  rate  of  ascent  differences 
which  were  considered  under  errors  1  and  2. 

f.  In  some  cases  when  a  reading  was  missed  at  one  installa¬ 
tion,  interpolated  readings  have  been  substituted.  Examination 
of  the  runs  indicates  that  this  can  be  done  with  sufficient  accu¬ 
racy  to  justify  the  practice. 

The  total  of  these  errors  i;a,  b,  c,  d,  e,  f  has  been  assumed  to 
average  3  per  cent  of  all  divergence  values. 


Table  6 

Individual  and  total  percentual  errors  for  various  divergence  values. 
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10 


Divergence  x  10"^  sec"* 
50  100 


150 


200 


1 

2 

3 

h 


23  5  3 

23  5  -3 

3  h  3 

3  3  3 


1  •  1 
1  1 

1±  ± 

3  3 


3.U 


Total 


33 


8.5 


6 


3.7 


31 


The  total  average  error  in  table  6  is  the  square  root  of  the  sum  of 
the  squares  of  the  individual  average  errors.  The  individual  errors  are 
assumed  to  be  independent  and  random.  This  total  average  error  is  graphed 
in  fig.  11  and  clearly  shows  the  dependence  on  the  magnitude  of  the  diver¬ 
gence.  The  mean  value  of  divergence  (for  the  first  twenty-seven  minutes) 
was  54.8  x  10"5  sec-*  and  this  value  is  about  8  per  cent  in  error.  Values 
smaller  than  20  x  10-5  sec-*  are  >20  per  cent  in  error.  The  vorticity 
computations  can  be  assumed  to  have  corresponding  errors.  The  mean  vorticity 
of  74.4  x  10-5  sec-*  is  in  error  by  8  per  cent  and  values  smaller  than 
30  x  10“5  sec-*  are >20  percent  in  error.  A  simulated  curve  for  vorticity 
has  been  included  in  fig.  11. 

8.  Conclusions 

a.  The  magnitudes  of  the  divergence  and  vorticity  measured  over  a 
small  area  of  about  £  square  mile  were  found  to  be  about  55  x  10“^  sec-* 
and  75  x  10-5  sec-*,  respectively. 

b.  It  is  felt  that  the  values  quoted  in  (a)  are  in  error  not  more 
than  10  per  cent. 

c.  Primary  and  secondary  levels  of  least  divergence  were  found  at 
6000  ft  and  13,000  ft.  In  [3]  similar  levels  were  found  at  8000  ft  and 
13,000  ft.  The  study  [3]  dealt  with  a  different  geographical  region,  a 
different  season  and  much  larger  scale  motion. 

d.  The  absolute  value  of  the  divergence  did  not  appear  to  vary  di- 
umally  (as  in  [3]),  but  there  is  some  evidence  that  it  varied  somewhat 
systematically  throughout  the  observational  period  of  about  six  weeks  (see 

fig.  4). 

e.  One  triple -balloon  run  showed  possible  evidence  of  a  very  narrow 
line  of  wind  direction  change  in  the  horizontal*  Postulating  the  existence 
of  such  a  line,  the  small  scale  wind  observations  can  be  shown  to  be  con¬ 
sistent  with  one  another  as  well  as  with  the  larger  scale  motion. 
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f.  It  follows  from  (e)  that  large  values  of  divergence  (+  10"^  sec“l) 
nay  be  found  locally  (over  a  £  square  mile  area).  However,  even  on  this 
small  scale,  such  large  values  were  found  to  be  the  exception  rather  than 
the  rule. 
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